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Abstract 

We calculate the deuteron anapole moment with the wave functions obtained from the 
Argonne v 18 nucleon-nucleon interaction model. The anapole moment operators are consid- 
ered at the leading order. To minimize the uncertainty due to a lack of current conservation, 
we calculate the matrix element of the anapole moment from the original definition. In virtue 
of accurate wave functions, we can obtain a more precise value of the deuteron anapole mo- 
ment which contains less uncertainty than the former works. We obtain a result reduced 
by more than 25% in the magnitude of the deuteron anapole moment. The reduction of 
individual nuclear contributions is much more important however, varying from a factor 2 
for the spin part to a factor 4 for the convection and associated two-body currents. 
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1 Introduction 



After Zel'dovich introduced the concept of the anapole moment (AM) [I], a first non-zero 
measurement was reported quite recently in the 133 Cs atom 0. Isolation of the effect is 
difficult because the AM has a small contribution compared to the leading order Z° exchange 
in atomic physics. However, the part that involves the nuclear spin is suppressed by a factor 
1 — 4 sin 6^ (~ 0.08). Thus in some cases, the AM which is a higher-order effect in electro- 
weak interactions can be comparable with the spin-dependent Z° contribution. Flambaum 
and Khriplovich showed that the AM of a heavy nucleus is proportional to A 2 ^ 3 ||, from 
which one can deduce the dominance of the AM over Z° exchange or radiative corrections 
for large A. Since then, calculations of the AM of heavy nuclei have been the object of the 
major interest in the theoretical works of the domain 0, [|, [|, [6|, [7[]. 

The dominant contribution to the AM of a light system, namely the deuteron, was 
pointed out by one of the author |J. The spin current which stems from the anomalous 
magnetic moment of the nucleon was expected to give a few times larger contribution than the 
convection current or exchange currents and it was verified in several papers || [10], [II], As 
a result, the contribution of the AM becomes similar to the radiative corrections. Khriplovich 
and Korkin (KK in short) calculated the AM of the deuteron analytically with the zero-range- 
approximated (ZRA) wave functions |JTTJ] . Their result of the spin-current term, which is the 
most dominant contribution to the AM of the deuteron, has exactly the same form as the 
one obtained from the framework of the effective field theory by Savage and Springer (SS in 
short) |l(J. The error of the result, which may be mainly from the simple wave functions, 
is estimated to be about 20%. However, in recent calculations of the asymmetry in the 
n + p — > d + 7, we showed that the ZRA (or the effective field theory [O]) result exceeds 



the ones obtained with a few phenomenological wave functions by more than 50% 0, [15 
With this observation, it may be possible that the error of the deuteron AM with ZRA wave 
functions can be larger than 20%. With the purpose to minimize theoretical uncertainties, 
we calculate the deuteron AM with the wave functions obtained from the Argonne v 18 
potentials p~6|| . 

The low momentum transfer which characterizes the process makes it possible to treat 
the problem with heavy-baryon-chiral-perturbation theory (HBChPT). With the counting 
rules of HBChPT, one can obtain the transition operators order-by-order in a well-defined 
way. The error of a calculation can be systematically estimated from the higher orders that 
are not taken into account in the calculation. Since well-ordered operators are evaluated 
with very accurate phenomenological wave functions, we expect that the uncertainties will 
be the least among the theoretical calculations. 
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2 The Anapole Moment 



The anapole moment of a system is obtained from the expansion of the vector-potential 

A,(x) = [ dx' r ^- (1) 



in a series in i 1 (x = |x|). The quantity j is the matrix element of the current density 
operator for given initial and final states 

j=<V/|j|V'i>- (2) 

The zeroth order term vanishes since there is no net current and the first-order term gives the 
vector potential of the magnetic dipole. The second order can be separated into a magnetic 
quadrupole term and the anapole term. After some algebra, one obtains the following form 
of the anapole term ^ 

aE i/ rfxxx ^ xx ^ 

with which the vector potential of the anapole moment reads 

A anapole (x) = — (-a V 2 ^ + a ■ V V-) . (4) 

The second term of A anapole can be removed by a suitable choice of gauge and the resultant 
anapole vector potential takes the form 

A ana p l e (x) = a 5 (3) (x). (5) 

If current conservation is satisfied, Eq. ([3D can also be written as 

a = — 7r J rfxi 2 j(x). (6) 

In many calculations of the nuclear AM, Eq. (B) was adopted as the working definition of 
the AM. However, as shown in [[/J, the contribution of a current term depends strongly on 
this definition while the total result may be sensitive on fulfilling the current- conservation 
constraint. The situation becomes more uncertain for the exchange currents or higher order 
terms. Thus in order to avoid possible large effect on the result due to lack of current 
conservation, we preferentially adopt in our calculations the original definition, Eq. We 
nevertheless stress that we include in our work a minimal set of two-body currents ensuring 
current conservation in relation with the parity-non-conserving (PNC) one-pion exchange 
interaction. This one (together with the PC part) turns out to be essential in ensuring the 
approximate equivalence of Eq. (6) with Eq. (3), as will be seen at the end of the paper. 



# 3 While we have the AM in fm 2 , the Seattle group favors a dimensionless definition which can be obtained 

N/ 



by multiplying ours by m 2 N /4ir 
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3 Operators and Wave Function 



i 



i i 

x x 

pair pion 

(a) (b) 

Figure 1: Diagrams representing one and two-body electromagnetic contributions. 

The diagrams considered are shown in Fig. 1. Leading one-body currents (Fig. l-(a)) 
are composed of spin and convection currents and PNC exchange contributions are in Fig. 



l-(b). The leading PNC vertex marked with x reads [17 



4^ = -^f Ni{?x*)'N. (7) 



The current density operators of each diagram are 



Jspin V 



J pair \X ) 6 



eE^ L V x x(^^)(x-r,,)), (8) 
eg^{ P ^ (3) (x-,)}, (9) 
(r x ■ r 2 - rfr|) y (r 12 ) £ a, t 5®(x - r,), (10) 



n h (1) 

jpion(x) = -e^^^-Ti-Tf^^.^-^-^^-^itoC^yoCrii.) (11) 
where r i2 = \t\ — r 2 |, r ix = |r» — x| and 

Vo(r) = 



Aur 

H l N is defined as 

with [1$ — 0.88 and \xy = 4.71. The PNC interaction of the proton and the neutron generates 
parity-odd components in the deuteron wave function. In the context of the meson-exchange 
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picture, the PNC interaction is mediated by 7r, p, to and heavier mesons. PNC components 
in the wave function can be obtained by solving the Schrodinger equation with the PNC 



potentials given in [pL7|1 - In a low energy process, it is believed that the pion exchange will 
dominate PNC interactions when its contribution is not forbiden by some selection rule. The 
pion-exchange PNC potential reads 



(12) 



where r = r p — r„, r = |r| and I = \{B P + a n ). The constants, g& and are given the values 
1.267 and 92.4 MeV respectively. For the deuteron, the PNC potential, Eq. ( |12|) gives rise 
to a parity-odd 3 Pi component. We write the parity admixed wave function as 



Airr 



u(r) + S 12 (v)^ff] Coo 



ih 



rNN\ 



I • r f (r) ( 
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(13) 



where Si 2 (f) = 3 (?i • f <r 2 • f — B\ • Oi and x an d C represent a spinor and an isospinor, 
respectively. The quantities u(r), w(r) and v(r) are obtained by solving the Schrodinger 
equation and one can calculate the matrix elements of the current density operators, Eqs. 
(8) - (11), with the obtained solutions. After that, the calculation of the anapole moment 
with Eq. (3) is straightforward. 

4 Results 



-"-spin 



elh 



(i) 



The anapole moment, Eq. (3), reads for each term 

~P L v\ -—— [ drrvir) (u(r) — \/2w(r)) 
V 6 m N J v ' 

-\ / drrv(r) (u(r) — \/~2w(r)) elh^ 

3 V 6 m N I ' 

V2irg A 



NNi 



*pair 



l pion 



9U 
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Numerical results are 



V2 

drry (r) ^u(r) + ^y^ w ( r )\ x 

u(r) (l - -m^rj - -^™( r ) + T^m^r 



elh 



(i) 

nNN- 



"spin 
a-conv 

^-pair 
^-pion 



-0.531 el h ( ^ NN , 



0.038 el h 



(i) 

:NNi 



-0.026 elh^fjf, 
0.027 el h^f N 



(14) 
(15) 



drr 2 y (r) lu(r) + -~w(r)) (u(r) - V2w(rfj elh^ ] NN , (16) 



(17) 

(18) 
(19) 
(20) 
(21) 
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Figure 2: The diagrams of the leading order nucleonic anapole moment. Solid, dashed and 
wavy lines represent the nucleon, pion and photons, respectively. 



where all the values are in the fm 2 unit. 

In the HBChPT, one can specify the magnitude of a diagram in terms of the power (u) 
of the momentum transfer (Q) by the rule 

v = 2L - 2(C - 1) - 1 + J2 v i ( 22 ) 

i 

where 

v i = d i + -£ + e i -2. (23) 

L is the number of loops, C is the number of disconnected lines, di, rii and are the numbers 
of derivatives, nucleon lines and external gauge fields at the vertex i, respectively. All the 
strong and electro-magnetic vertices satisfy 

Vi > (24) 

but the PNC vertex, Eq. (7) has 2/j = —1. One can easily verify that the diagrams in Fig. 
1 have v = —2. Nucleon anapole diagrams shown in Fig. 2 also have this order. In order to 
be consistent with the order counting of the effective field theory, the nucleonic terms also 
should be included. 



This term was already calculated in several works |L8|, [19], |2(| |2IJ . Since the results are 
consistent in both magnitude and sign, we do not repeat its calculation but just adopt the 
result here. The leading order nucleonic anapole, taking into account the deuteron D-state 
probability Pd, reads 



6V2/ 7r m 7r 2 



a N = (l--P D )e!K 

r 

-0.417 e I (25) 
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Then the anapole moment of the deuteron at the leading order is 



a d a spin <^conv a pair a pion ~l~ a N 



-(0.531 - 0.038 + 0.026 - 0.027 + 0.417) elh 



(i) 

'ttNN 



0.909 e I h£l N . (26) 



5 Discussions 



As expected, the spin term is the most dominant one among the contributions from the 
currents. The contributions from the convection, pair and pion terms are much smaller 
while the two last ones strongly cancel. The sum of these three contributions is only about 
7% of the spin term. However the nucleonic term is close to the spin term, which is quite 
contradictory to the former results of the deuteron anapole moment. 

In order to illustrate this feature clearly, let us compare our result with the ones by SS 
Til an d KK [11]. Firstly, SS's result is 



a cF — a spin + a PE + a N 

= -(1.03 - 0.18 + 0.46) e Ih$ fN = -1.31 e I h^ NN , (27) 

where the pion-exchange term (PE) corresponds to the sum of the convection, pair and pion 
terms in our calculation. We have verified in a separate calculation that SS's result can be 
obtained with the ZRA wave function P2|. Secondly, KK's result reads 



a KK _ „KK , „KK , _KK 
d d — ^spin d orb "r d N 

= -(1.03 -0.07 + 0.21) e!h^ NN = -1.17elh^ NN . (28) 

Since KK also used the ZRA wave function, their spin term is equal to SS's one. However 
there are substantial differences in the remaining terms. In the following, we successively 
discuss these different contributions. 

The biggest quantitative difference of our result with the above ones comes from the spin 
term. The Avl8 result is roughly half of the ZRA value. In ||11|| , KK argued that the uncer- 
tainty of the ZRA wave function is about 20%. Taking this uncertainty into consideration, 
the spin term can be reduced to —0.82 el h^ N which is still larger than our result by about 
55%. The discrepancy of our result with the KK (as well as SS) estimate is understood by 
noting that our calculation incorporates the effect of a short-range repulsion in S-states as 
well as the known repulsive character of the NN interaction in the 3 Pi channel. We also note 
that the ZRA wave function contains only a central component, u(r), while our calculation 
includes a tensor one, w(r) too. Contrary to the PNC asymmetry in n-p radiative capture 
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T4| , p~5|1 , this component produces a partly destructive interference. The contribution of each 



part is 

= -0.724 e I (29) 
aj, = 0.193 e I (30) 

In some cases, the contribution of the tensor part in the wave function is assumed to be small 
and its contribution is neglected. However in the deuteron anapole moment, this contribution 
is non-negligible and its effect on the magnitude of the anapole moment is substantial. 

We now turn to the second contribution in Eqs. (|27| , |28|) , which, among other contribu- 
tions, involves the convection current one. There too, a large suppression of our results is 
observed but the effect is more drastic for the PE term in SS's result than for the orbital one 
in KK's result. In SS's result, PE term's contribution is about 15% of the total magnitude, 
but in our result it is only 4% of the total value. Then, the question arises why the PE term 
in SS does not coincide with the orbital term in KK. 

It can be shown that the anapole moment of the convection term with the definition Eq. 
(3) is equivalent to the orbital term 



TiC 



-I- 



l 2 



(31) 



12m7v 

where 1 = r x p is the angular momentum operator in the center of mass frame. The global 
factor of the sum of the spin and the convection term, (/j, v — 1/3) always appears regardless 
of wave functions. This explains the relative ratio of the contributions due to the spin and 
convection currents in KK's work as well as in ours. 

As mentioned above, KK's orbital term accounts for only the convection term. In KK's 
argument, they made use of the fact that 

Jpair(r) + Jpion(r) oc r. (32) 

With the definition Eq. (3), it may be that this term gives a zero contribution since 

aocrx(rxj). (33) 

We would like to point out that caution is required in the above reasoning. Firstly, the 
matrix elements that are relevant in the calculation of the anapole moment involve the 
current density operator and not the current operator. It can be easily shown that the 
current operator which is obtained by integrating the current density operator with respect 
to the field point x does satisfy Eq. fl32l). However the anapole moment should be obtained 
from the double vector product of the current density operator and field point x, whose 
result may not be like Eq. (|32|) in general. Secondly, a rough substitution x — > r in Eq. (3) 
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gives Eq. fl3"3"|). This substitution, even if one is working in the center of mass frame and 
thus center of mass coordinate is discarded, should be carefully derived from the evaluation 
of the matrix elements with the two-body wave function. A crude transformation of the 
coordinates may give wrong results. As our result shows, a rigorous derivation with the 
current density operator gives a non-zero contribution of the two-body currents. We also 
checked that this derivation provides coincidence with SS's result |[22|| . Curiously, the extra 
contributions that the SS's work accounts for do not seem to show up in our results. It turns 
out that the term proportional to l/m^ in Eq. (0), a priori favored, is suppressed with the 
ZRA value (numerical factor 0.07 instead of 0.18). It is reminded that in the zero-pion-mass 
limit, this term is cancelled by another one arising from the nucleonic term [K|. The other 
term (Eq. (^) and part of Eq. fl2T|)), which is less singular in configuration space and has 
an opposite sign, is essentially unchanged. As a result, the sum of the two contributions 
almost vanishes in our work. 

Concerning the nucleonic term (third one in Eqs. ([27], |28|)), KK adopt 



The first term in the parenthesis coincides with the one we retained and the second term 
stems from the correction to the first term. In the context of the counting rule of 

the effective field theory, this l/m^ term is classified in the higher order corrections to the 
leading term. In order to be consistent with the strategy of the effective field theory, the 
1/mjv corrections should be taken into account consistently, i.e. their correction should be 
calculated not only for the nucleonic term but also for current terms or any other types that 
have the same order. In that sense, KK's result contains the l/m^ correction partially. How- 
ever, it is interesting to notice that the magnitude of the l/m^ correction to the nucleonic 
term amounts to a half of its leading value. This indicates that the higher order corrections 
can play a critical role in the magnitude of the deuteron anapole moment. Their calculation 
should be done in the future. 

In a separate work, we have observed that, while the spin current satisfies gauge invari- 
ance by itself, convection, pair and pion currents are not gauge invariant solely and gauge 
invariance is restored when the three terms are summed up in a minimal case [Q. In the 
general case, one should also include two-body currents related to the description of the NN 
strong interaction. As a way to investigate how much gauge invariance is broken, we compare 
the anapole moments calculated with Eq. (3) and Eq. (6). Equation (3) is defined from the 
definition of the anapole vector potential and Eq. (6) is derived from Eq. (3) by imposing 
current conservation. Calculation of the anapole moment from Eq. (6) is straightforward 
too. The results are 




(34) 




0.051 e I K, 



■ttNNi 



(35) 
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(36) 
(37) 



where the superscript CC implies the result with the assumption of current conservation. 
As expected, individual contributions differ from those in Eqs. (p^-|21|), but it is also noticed 
that the sum differs, 0.020 e I h!$f N , instead of 0.039 eIh^ N . The first candidate to explain 
the discrepancy is the two-body current associated to the strong one-pion exchange inter- 
action. The corresponding contribution, 0.018 e I /4aw> which should be added to the sum 



of contributions of Eqs. (|35| - |37|) , removes most of the difference. This one-pion exchange 
interaction can also generate self-gauge invariant contributions that will affect equally the 
sum of contributions in Eqs. ( ]T9| - |2~T|) and Eqs. (|35| - p7D . Due to cancellations of pair and pion 



contributions, this common contribution is found to be small (— 0.0007 el h^ NN j 

Concluding, we calculated the anapole moment of the deuteron with the wave functions 
obtained from the Avl8 potential. Its magnitude is reduced by more than 25% from the 
previous ZRA results but is still comparable to the radiative corrections. If the contribution 
of the nucleon anapole moment is put aside, the effect is much larger however, ranging 
from a factor 2 for the spin contribution to a factor 4 for the contribution of the convection 
current and associated two-body currents. We observed that gauge invariance for these last 
contributions, hopefully smaller, is a severe constraint and its fulfillment is a non-trivial 
problem. In the analysis of other's work, we noticed that the corrections of higher orders 
can modify significantly their estimate. The calculation of gauge-invariant higher-order 
corrections will be the future challenge. 
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